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Abstract 

In the confining phase of any gauge system the mean squared width of the colour 
flux tube joining a pair of quarks should grow logarithmically as a function of their 
distance, according to the effective string description of its infrared properties. New 
data on 31? 7L2 gauge theory, combined with high precision data on the interface 
physics of the 31? Ising model fit nicely this behaviour over a range of more than two 
orders of magnitude. 
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1 Introduction 



The mechanism of confinement for abelian and non-abehan gauge theories in three or four 
space-time dimensions seems well described by the old conjecture [Q] that the vacuum 
behaves as a dual superconductor. The dual Meissner effect squeezes the colour field 
generated by a pair of quark sources into a thin flux tube (the dual version of the Abrikosov 
vortex), generating a static potential proportional to the interquark distance. Numerical 
studies of the distribution of the colour flux around the quark sources ||2|, |3|, Q seem to 
support this picture. 

According to another old conjecture ^, this flux tube can vibrate as a free 
string. There are many interesting and impressive quantum effects associated to these 
string vibrations. One of them seems to conflict with a property of the dual Abrikosov 
vortices: the intrinsic thickness of these is a physical quantity which is independent of the 
interquark distance and is fixed by the parameters of the Ginzburg-Landau formulation 
of the superconductivity, while in the string picture the width of the flux tube grows 
logarithmically as a function of the interquark distance [Q. 

In this paper we will show compelling numerical evidence of this logarithmic growth. 
However this does not mean that the picture of the dual Abrikosov vortex does not work. 
Simply, the way one defines the mean squared width of the flux tube in the lattice sim- 
ulation also takes into account the quantum effect of the vibration modes of the dual 
Abrikosov vortex, yielding an effective width which is larger than the intrinsic thickness 
of a non-vibrating vortex. 

It is worthwhile to stress that the dual superconductor nature of the confining vacuum 
and the string-like behaviour of the colour flux tube are conceptually two independent 
properties of the theory. Indeed in the strong coupling regime, before the roughening 
transition, i.e. in the region where the strong coupling expansion of the string tension a 
converges, the dual Meissner effect is already operating and the colour flux tube connecting 
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pairs of quark sources has a string-like form, however cannot vibrate. Hence in this strong 
coupling phase the dual Abrikosov vortices, like those of the ordinary superconducting 
media, do not exhibit any kind of quantum string effects. 

At the roughening point, the colour flux tube of whatever 3-0 or 4Z) gauge theory 
undergoes a transition towards a rough phase (which is the one connected to the continuum 
limit), where it can vibrate freely and the energy stored in its vibration modes is no longer 
negligible. The collective degrees of freedom describing such a string-like motion can be 
described by a 2D field theory known as effective string theory. 

It is widely believed that the roughening transition is described by the universality 
class of Kosterlitz and Thouless Q. Accordingly, the renormalization group equations 
show that in the rough phase the effective string theory describing the dynamics of the 
flux tube flows at large scales towards a massless free field theory. Thus, for large enough 
interquark separations, it is not necessary to know explicitly the specific form of the 
effective string action describing the behaviour of the flux tube, but only its infrared 
massless limit. This simple, general fact has some important implications because the 
limit action is critical in the whole rough phase. As a consequence, it produces universal 
(i.e. gauge group independent) finite size effects. Moreover, because this theory is also 
infrared divergent, such finite size effects are expected to be rather strong. 

The effects we shall describe require highly accurate numerical data that at present 
can hardly be reached in 4D SU{2) and SU{3) gauge models. 

Luckily, these effects can be checked accurately in the 3D TL^ gauge model which is 
the simplest prototype of gauge theory. 

2 Width of the Colour Flux Tube 

The logarithmic behaviour of the mean squared width 'uP'{R) of the colour flux tube as 
a function of the interquark distance R has been predicted many years ago by Liischer, 
Miinster and Weisz in the framework of the effective string picture of gauge systems. 
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We shall rederive and refine this universal law by directly using some exact results on 
the two dimensional free gaussian model in a finite box and compare it with Monte Carlo 
simulations on the 3D 2L2 gauge model. 

We shall see that the free gaussian model describes quite accurately this phenomenon. 

2.1 Results from the Gaussian Limit 

In order to fix notations, let us recall that the simplest way to implement the effective 
string picture of the colour flux tube is to assume that the vacuum expectation value of a 
large Wilson loop W{C) can be represented by the functional integral 

{W{C))= Jl[[Dh]e-I''"^^^''^^ . (1) 

1=1 

For simplicity C is taken planar and a = 1, 2 are the coordinates of the planar domain 
V bounded by C. are the transverse coordinates of the string describing the 

configurations of the flux tube in its space-time evolution. They are the fields entering in 
the effective string action S = /d^^£(/ii), which is largely unknown, but is believed to 
flow to the free gaussian model for very large loops C, as explained in the introduction. 
More precisely we have 




where a is the string tension and A = Jj^ d?^ is the area of the planar domain. Eq. (Q) can 
be considered as the first two terms of an expansion in the adimensional parameter 1/aA. 
Recently, the first non-gaussian correction of such an expansion has also been studied 

The mean squared width of the flux tube generated by this Wilson loop is deflned 
as the sum of the mean square deviations of the transverse coordinates hi{^^,^'^) of the 
underlying string, i.e. 

wHe,e) = E {{Me,e)-hf'^y) , (s) 

i=l 
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where hf^ is the transverse coordinate of the center of mass of the flux tube. In terms 
of the Green functions Gi we get 

D~2 



(4) 

where e is a UV cut-off, a is the string tension and the Green function is defined as 

G^i^,e)=^ {h^{Oh^{e)) ■ (5) 

The vacuum expectation value is taken with respect to the string action S. In the infrared 
limit the Green function fulfills the free field equation 

-AG(e,0 = <^^'^(e-0 , (6) 

where for simplicity the index i of the transverse direction has been omitted. Instead of 
considering the usual expansions of G in terms of the eigenfunctions of A it is possible 
(and much more useful) to write it in a more compact form. 

Indeed the problem of finding the Green function for an arbitrary, simply connected 
region V can be solved in a closed form once a conformal mapping z = + i — > £ 
of V onto the unit circle \i\ = 1 is found which maps z' = ^'^ + i^'^ into the origin 
i = 0. Denoting by iz'i^) the analytic function providing us with such a mapping, it 
is immediate to verify that the real function /^/(z,!?) = log|^2'(2)| is harmonic in the 
punctured set D \ {z'} (hence A/ = for z ^ z'), vanishes at the boundary dD and 
diverges logarithmically as log\z — z'\ for z ^ z' . It follows that the Green function is 
given by 

Gv{z,z') = -^fAz,V) . (7) 

Some explicit examples can be found in Tab. I and will be discussed later. 

Inserting such expressions in Eq.(^), we can evaluate analytically the shape of the flux 
tube as a function of the coordinate z inside the domain D. One can observe a very steep 
increase of w'^{z) at the boundary of V followed by a plateau covering approximately the 
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of the domain, where w'^{z) grows very slowly an reaches its maximum in the center 
Zo of the domain. For instance, in the case of an infinite strip of width R we get 



2TTaw (r) 



log 



sin(2r + ^) 
|sin(f )| 



(8) 



where r = ^x, < x < R, and the UV cut-off has been reabsorbed in the length scale 
a measure of the total width of the flux tube it is often used in the literature 
the integrated quantity 

w' = jj^d'^w\e,e) ■ (9) 

However in the applications to the gauge models this quantity may be affected by non- 
universal corrections due to the self-interactions of the quark line at the boundary. For 
this reason we prefer to use instead the value of the mean squared width measured at the 
symmetry point of the domain: 



(10) 



For instance, in the Wilson loop of size R x T we have obviously = "f ?o ~ 



Tab. I 



V 



b.c. 



Green function 



Disk fixed G d{z, z') = —-^log 

Square fixed Gr{z, z') = —-^log 

Strip fixed G s{z, z') = -^\og 

Torus periodic Q{z — z' 



R 



(t{z-z')c7{z + z') 

a{z-z')cj{z+z') 

s\n{-K{z-z')/2R) 
sin(7r(z+2')/2iJ) 



0.8105 
-0.0139 

0.8587 




The closed form of the Green functions for some simple domains is given together with 
the difference between the mean squared width as defined in Eq.s ({^[7^ and that relative 

to a square torus. 
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For our purposes, the relevant property of the conformal mapping ^ is that one can per- 
form an arbitrary scale transformation z — > Az without destroying the conformal character 
of £. More precisely we can write 

U{z,V) = fA,,{Az,VA) , (11) 

where Pa denotes the scaled domain. A direct consequence is that it is always possible 
to fix the area of the scaled domain Da to an arbitrary value, say 1, without changing 
the Green function. It follows that the integration of the finite part G{z, z') in Eq. 
cannot depend on the size of the domain D but only on its shape. On the contrary the 
UV divergent part gives a contribution which grows logarithmically with the size of the 
domain. This can be simply understood as follows: if ii is a typical linear dimension of 
the domain D, the scaling property described above implies that the cut-off should appear 
always in the ratio e/R. On the other hand, for e « R the Green function behaves like 
— 2^1og(e), yielding the logarithmic law 

wl = ^log{R/R,) , (12) 

where the UV cut-off has been absorbed in the definition of the scale i?c. 

Note that the above equation is meaningful only for R > R^. This inequality has a 
simple physical meaning: the dynamics of the flux tube is described by a truly free field 
theory only at large distances; the cut-off e sets up the ultraviolet scale Rc below which 
either the free-field approximation breaks down or the internal degrees of freedom of the 
flux tube become important. Thus Rc gives us a very rough estimate of the intrinsic 
thickness of the flux tube. 

The absolute value of Rc is not directly calculable, being proportional to the UV cut- 
off. We can however evaluate the ratios among the Re's of different domains making the 
very mild assumption that the UV cut-off e cannot depend on the size and the shape 
of the T> domain (which are on the contrary infrared properties). For later convenience, 
we shall use as reference scale Rp associated to a squared box with periodic boundary 
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conditions {i.e. a torus). Also in this case the the Green function Q can be expressed in a 



closed form. More generally, for a torus of sides RxT we have (see for instance Ref . [ 11 1 ) : 
G(z) = -^ln\a{z)\ + ^^^^M^') + ^{Qmzy , (13) 
where ct{z) is the Weierstrass sigma function, defined through the infinite product 



a(z) = 2 TT ( 1 - - ) , uj = mR + inT . (14) 



The set {uj} is known as the period lattice and E2{it) is the first Eisenstein series with 
t = T/R. For more details see also the Appendix. 

Denoting with Wp{R) = i^\og{R/ Rp) the mean squared width of the flux tube 
generated by taking T> to be the square box of reference, it follows that the difference 
Aw'^{R) = w1{R) — Wp{R) does not depend on R but is a calculable function of the shape 
of P. 

In Tab. I one can find the closed form of the Green function and the corresponding 
Aw'^ of some simple domains like a disk of radius R, a square of side R, an infinite strip 
of width R and the reference torus of size R x R {R = T). It is important to notice 
that the Weierstrass sigma function entering the definition of the Green function of the 
square of side R corresponds to a torus of double size. For a rectangle of size RxT the 
Green function has the same formal expression as the square in terms of cr{z), which now 
is defined on a torus of size 2R x 2T. 

The detailed calculation of w'^ for the torus (which is interesting in the physics of fluid 
interfaces) and for the rectangular Wilson loop can be found in the Appendix. The result 
is drawn in Fig. 1, where the behaviour of Aw'^ for a rectangle as a function of the aspect 
ratio t = T/R (dashed line) is compared with that of a rectangular torus with the same 
shape t (continuous line) and with the value for the infinite strip (dotted line). Note that 
the almost linear increase of Aw'^ for the torus agrees with the fact that for large t the 
theory becomes almost one-dimensional and describes the interface of a 2D Ising model, 
where the mean squared width is known to grow linearly with the size. 



This behaviour of IS.w^ calculated in the free gaussian model as a function of the shape 
and of the boundary conditions (periodic and fixed) can be easily checked in the 2D XY 
model (in Fig. 1, open circles and triangles respectively) which is believed to belong to the 
Kosterlitz-Thouless universality class and certainly flows to the free gaussian model at 
large scales. The other numerical data reported in Fig. 1 (black dots) are data extracted 
from the SL* ^2 gauge model simulations which will be discussed in the next section. 

If one is able to evaluate, besides the effective squared width, also the higher moments 
^(2n)^ defined by 

^^'''^ = \i: j d'i{{h,{i\e)-h<t'')"'') , (15) 

one can retrieve the three-dimensional distribution of the colour field strength density 
p{x, y, z) inside the fiux tube. Indeed, choosing for simplicity the domain T> sitting in the 
plane z = 0, we can implicitly define such a density through the following equations 

^(2n) ^ Jpz^'^dx dy dz ^ ^^^^ 
f pdxdy dz 

As a consequence, the Fourier transform of p can be written in terms of these moments as 

00 2n 

p(p) = E(-ir7f3T^^'"^ ' (17) 

n=0 ^ ' 

where the integral of the density has been normalized to 1. In the free gaussian model the 
above calculation can be performed explicitly using the Wick theorem, yielding simply 

tt;(2n) ^ (2n - 1)!! (w;2)" . (18) 

Inserting this result into Eq. (|l7|) one finds that the density has a gaussian distribution 

Pg{z) = aex.-p{-z^/b^) . (19) 
2.2 Simulations on TL2 Gauge Model 

Though the logarithmic growth of the squared width of the flux tube is the most important 
and model-independent quantum effect predicted by the effective string description, it has 
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not been observed until now, because it is an infrared phenomenon that can be seen only 
in very large Wilson loops, which are at the limit of the sizes reached by present numerical 
simulations on SU (2) and SU (3) 4D gauge theories. However, it is nowadays possible to 
overcome this problem in the 3D gauge model by exploiting some special features of 
this model. 

The 3D 7L2 gauge model on a cubic lattice is defined through the partition function 

Zgauge{P)= E ^xp , (20) 

where the action 5 is a sum over all the plaquettes of the cubic lattice 

S = -'^aa , aa = aij^ai^ai^ai^ . (21) 
□ 

This model can be translated into the usual 3D Ising model by the usual Kramers-Wannier 
duality transformation 

ZgaugeiP) ^spiniP) (22) 

^ = -^log[tanh(/3)] , (23) 

where Zgpin is the partition function of the Ising model on the dual lattice 

ZspiniP) = E eM-PHis)) , (24) 

Si = ±l 

with 

H{s) = -^J{ij)SiSj (25) 

where i and j denote nodes of the dual lattice and the sum is extended to the links (ij) 
connecting the nearest-neighbour sites. For the moment the couplings J^j^^ are all chosen 
equal to +1 . 

Using the duality transformation it is possible to build up a one-to-one mapping of 
physical observables of the gauge system onto the corresponding spin quantities. For 
instance, the vacuum expectation value of a Wilson loop W{C) can be expressed in terms 
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of spin variables as follows. First, choose an arbitrary surface S bounded by C: 9S = C\ 
then "frustrate" the links intersecting S, i.e. take J^j^^ = —1 whenever {ij) n S 7^ 0. Let 
us denote with H'{s) the Ising Hamiltonian with this choice of couplings: the new Ising 
partition function Z'gp^^{j3) = Yl,Si=±i ^^P ("/^-^^'(•s)) describes a vacuum modified by the 
Wilson loop W{C), which we shall call the W-vacuum. A well-known consequence of 
duality is that 

Z' 

{W{C = dT.)) gauge = ^ = ( H exp(-2/3s,S, , (26) 

where the product is over all the dual links intersecting S. 

One of the main difficulties in dealing with numerical simulations is that near a critical 
point the successive configurations generated by a Monte Carlo procedure are strongly 
correlated. Then it is hard to get a reliable set of thermalized configurations where the 
measurements can be done with sufficient accuracy. This phenomenon is known as critical 
slowing down. 

A great advantage of mapping a gauge observable onto a spin observable is that a 



non-local cluster updating algorithm |12] can be used in place of the usual Metropolis or 
heat-bath methods. This kind of algorithm substantially reduces critical slowing down. 

This allows us to probe the structure of the flux tube with a high accuracy. There is 
no analogous procedure for the other gauge theories. 

The procedure is the following. A Wilson loop W{C) is realized in the spin lattice by 
frustrating all the dual links intersecting a given surface S bounded by C. These frustrated 
links modify the vacuum state so that the expectation value {P)w of the plaquette, or 
better its spin counterpart described in Eq. (]26|), becomes a function of its relative position 
with respect to W{C) and is related to the expectation value in the ordinary vacuum by 

{P)w = {W{C)P)/{W{C)) . (27) 

The difference between the expectation value of the plaquette in the vacuum modified 
by the presence of W{C) and in the ordinary vacuum can be considered as a measure of 
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the colour flux density of the flux tube. Choosing for instance as a probe a plaquette Py 
parallel to the plane of the Wilson loop we can take 

p\^{x,y,z) = {P\\)w-{P) . (28) 

Other orientations of the plaquette give approximately the same distribution (this is not 
the case for other systems 

It is important to stress that within our procedure we do not need to evaluate directly 
the correlation function {W{C)P) between the plaquette and the Wilson loop, which 
is very difficult in ordinary lattice gauge model simulations and must be handled with 
special techniques like the cooling 1^] or the smearing methods Q. Actually we only have 
to evaluate, according to Eq. (|28|), the mean expectation value of the plaquette with two 
different actions: the ordinary action and the one modified by the Wilson loop. This allows 
us to estimate p with high accuracy even with relatively modest computational means. A 
preliminary account of these numerical simulations has been reported in Ref. [T^ ]. 



Tab. II 





lattice size 


R 


T 


# data 


y 

Wo 


0.7516 


24x24x48 


12 


12 


12500 


19.0(1.1) 


0.7460 


22x22x40 


11 


11 


35000 


12.3(5) 


0.7460 


30x30x48 


15 


15 


7000 


15.6(8) 


0.7460 


60x60x48 


30 


30 


5000 


20.3(4) 


0.7460 


120x120x32 


60 


60 


1000 


24.8(6) 


0.6543 


64x64x32 


32 


32 


7400 


2.57(2) 


0.6543 


124x124x24 


32 


64 


2000 


2.79(5) 


0.6543 


96x192x32 


32 


96 


1700 


2.75(5) 


0.6543 


124x124x20 


64 


64 


1000 


2.95(3) 



The first two columns show the f3-values and the corresponding size of the lattices used 
in the present study, in the notation Lx x Ly x where R and T point in the x and y 
directions respectively. In the last four columns the size of the Wilson loops, the number 
of data and the corresponding estimates of the mean squared width of the flux tube are 

reported. 
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We would also like to point out that the plaquette expectation values entering into 
Eq. (p8|), namely {P\\)w (with the W-vacuum) and (P) (with the ordinary vacuum), 
are evaluated in distinct runs, so that the corresponding measurements are completely 
uncorrelated. 

We analyzed three distinct values of the coupling (3 and different sizes and shapes of 
square and rectangular Wilson loops, as can be seen from Tab. II. 

For each run with a W-vacuum, the Wilson loop is located in the zq plane. Cross- 
correlations in the measurements of the plaquette expectation values defined at different 
distances Iz — zqI from the Wilson loop have to be carefully considered in this procedure. 
To keep them under control we did five update sweeps of the lattice between two measure- 
ments, and we measured only plaquettes in a fixed z — zq slice each time. This means, for 
instance, that two measurements of the same observable were separated by 5 x sweeps, 
and there were at least 5 sweeps between measurements of different observables. From 
Tab. II one can see that each run contained between 100 000 and 240 000 sweeps. 

Moreover, if one takes a plane parallel to the Wilson loop and looks at the shape of 
p as a function of the x-y coordinates on this plane, one can see that p is approximately 
constant, except in a region extending over a few lattice spacings from the border of the 



Wilson loop [14|. Thus, in order to have a sufficiently stable evaluation of the width at 
the center of the Wilson loop, we averaged the plaquette expectation values {P\\)w over 
the x-y directions, excluding a border along the Wilson loop perimeter of size four lattice 
spacings across. 

A typical shape of the distribution p measured on the symmetry axis of the loop 
as a function of the distance z — zq from the loop plane is drawn in Fig. 2. The dotted 
straight line corresponds to the plaquette expectation value in the ordinary vacuum, {P) = 
0.8945(1) for (3 = 0.7460 (while we obtained (P) = 0.70852(2) at /? = 0.6543 and (P) = 
0.9108(2) at /3 = 0.7516): it does not show appreciable finite size effects for the lattices 
we have considered. The other two curves will be explained in the last section. 
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The mean squared width of the flux tube can be defined as 

2 / z^pi^\dxdydz ^^^^ 
° / p|| dx dy dz 

The integration region in the horizontal directions is taken over the central plateau where 
p is approximately constant as a function of x and y. The evaluation of the integral was 
performed numerically and all the errors were evaluated with the binning procedure and 
the jackknife method. 

It is important to stress that the numerical integrations are not reliable if the asymp- 
totic values of p and (P) do not agree within errors. To this end it is crucial to take large 
lattice sizes compared to the Wilson loops (see Tab. II). 

An important feature of Eq. (|l2|) is that it can be written in a universal, /3-independent 
form by expressing all the dimensional quantities Wq, R and Rc in units of y/a. Accordingly, 
we report in Fig. 3 the data of Tab. II and also very accurate data |jl^ for the mean squared 
width of fluid interfaces in the Ising model (black dots) for one of the couplings analyzed 
in our approach (/3 = 0.6543 of the gauge model). Note indeed that, according to Tab. I, 
the squared width of the flux tube with periodic boundary conditions (which is the case 
of the fluid interfaces) is almost the same of that with fixed b.c. 

The straight line represents a fit of all the data to the two-parameter formula 



= a 



log ^/^R + b . (30) 



We used the values a = 0.0107(1) at (3 = 0.7516 and a = 0.0189(1) at (3 = 0.7460 



according to Ref. |11|. Far from the scaling region, the string tension should be replaced 
by the stiffness k. This is the case for (3 = 0.6543 where we used k = 0.234(9) according 
to Ref. 0. One gets b = 0.173 ± 0.009 and a = 0.150 ± 0.005 in good agreement with 
the theoretical value ath of the parameter a, fixed by Eq. (|l^) to be a^^ = ^ — 0.159155 . 

Within high statistical accuracy these data clearly support a logarithmic widening of 
the flux tube in a range of quark separation R over more than two orders of magnitude, 
starting at about ^/aRmin — 1-2. 
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The physical adimensional quantity \faKc is expressed in terms of a and h as \faRc = 
Q-b/a^ Using the fitted values of the parameters we get y/aRc = 0.32 it 0.02. 

In order to compare these results with analogous data for other gauge groups, we may, 
with an abuse of language, express y/a in the physical units of the QCD string . Then the 
infrared logarithmic behaviour observed in Fig. 3 starts at least at Rmin — 0.6 fm while 
the maximal probed elongation of the flux tube corresponds to more than 100 fm. The 
data on 4D SU{2) flux tubes ^ cover now a distance up to 2 fm, but are still affected 
by strong systematic errors and are compatible also with a constant width for distances 
larger than 1 fm. 

3 Discussion 

The observed logarithmic widening of the flux tube exhibits, for the first time, a macro- 
scopic effect due to the string-like motion of the colour flux tube. At present it is visible 
only in the 3D TLi gauge system because only there it is possible to reach, using scaling and 
the duality property of this system, large separations of the quark sources. Nevertheless 
we expect a similar behaviour in all 3D and 4D gauge systems, because the appearance 
of string-like vibration modes of the colour flux tube is a universal property of the rough 
phase of whatever gauge model. 

The agreement of the logarithmic growth of the mean squared width of the flux tube 
observed in TL2 gauge model with the predictions of the gaussian model seems rather im- 
pressive: indeed the universal slope of the logarithmic growth has been checked in a range 
(more than two order of magnitude) which is unusually wide in numerical simulations. 

Admittedly, this prediction is not specific to the gaussian model: such a model is a 
member of a one-parameter family of conformal field theories, namely those with confor- 
mal anomaly c = 1. These theories have the same local behaviour (hence produce the 
logarithmic growth of the mean squared width) but have different topological excitations 
described by winding modes around a circle of radius r. This is the parameter labelling 
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the different members of the family and r — > oo corresponds to the gaussian limit. 

Also the agreement with some predictions which are more specific of the gaussian 
model is rather good. In particular, the scale Rc of the logarithmic growth depends also 
on the shape t = T/R of the rectangular Wilson loop. As a consequence, the difference 
Aw'^iR) = w1{R)-wl{R) between the mean squared width of the flux tube generated by 
the rectangular Wilson loop and that of the reference torus does not depend on R but is 
a calculable function of the shape of t as shown in Fig. 1. For fixed boundary conditions 
(dashed line) it reaches the asymptotic value of the infinite strip for very large values of t. 
The data extracted from numerical simulations in 'Sj2 gauge model (black dots) seem to 
agree with this behaviour, even if the errors are too large to reach a definite conclusion. 
The fact that the data seem a bit lower than the expected values is probably due to a (non 
universal) finite size effect; indeed in the 2D XY model with fixed b.c. such an effect can 
be easily observed. 

Another specific prediction of the gaussian model is that the scales Rc for fixed and 
periodic boundary conditions are almost the same. This fact is nicely confirmed with high 
accuracy by Fig. 3. 

The only feature of the gaussian model which seems to slightly deviate form the ob- 
served data concerns the distribution of the flux density as a function of the distance 
from the Wilson loop. It should be, according to Eq. (|l^), a gaussian Pg{z). Although pg 
(dashed line in Fig. 2) reproduces the raw features of the data, the fit is not very good. A 
much better fit is obtained by putting p = Pg/{1 — "fPg), where 7 is another free parameter 
(continuous line). This is suggested by the assumption that the flux tube does actually 
self-interact and the interaction should be a power series in p. 

In conclusion, the above results seem to indicate that the free gaussian model is a good 
approximation of the string-like behaviour of the flux tube at least in the infrared limit. 
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Appendix 

The Green function of the 2D free massless theory on a rectangle of size R x T with 
periodic boundary conditions on both sides, as given in Eq. ([l3|), can be expanded as |17] 



with 



0* = ^§i^^2'('') . t = T/R (A.2) 



where C is the Riemann zeta function and the Eisenstein series have the following expansion 
i?2fc(T) = l + (-l)^-f E^2fc-iHg" , (? = e2-- . (A.3) 

n=l 

(7fc (n) denotes the sum of A;th powers of the positive divisors of n and Bk are the Bernoulli 
numbers. 

Since Q(z,z') is translationally invariant under z ^ z + a; z' — > z' + a, it is only a 
function of the difference z — z'; then the mean squared width as defined in Eq.(^) w'^ is a 
constant with respect to z that we denote with Wpf.j.{R). Performing the integrations we 
get at once 

2TTaw1^j.{R) = log [r^/1 + t^/2e^ + arctant + | arctan 7 - | - f| 

- T.T=i ^'£luT^^ ((1 + t^)^^tB, sin {{2k + 2) arctan t) . (A.4) 

This expansion converges only for VT+~i^ < 2 because there is a pole sX z = 2i-K in the 
complex plane of the variable z = ■n{R + iT)/R. However it is possible to perform an 
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analytic continuation which allows us to evaluate nP' for any value of t, by resumming the 
potentially divergent part of the series. Indeed exploiting the well-known expansion 



oo rj 

1 2 f^} ' 2k] 



k=l 



we can easily derive the result 



^tt^ ' {2k + 2)1 2k 



B. z2fc+2 oo^ ^-nz 



E 

n=l 



1 z\ tt' 
- + ylog. + -.-- 



So z 
z^ 



3 

12 



(A.5) 



(A.6) 



which is useful for large \z\. Its insertion in Eq. (A^) gives the wanted resummation, once 
one realizes that E2kiit) = 1 + 0(e~2'^*). In particular, for large values of t one gets the 
asymptotic expansion 



vr R 
2Trawl^^{R) ^ -t + log— , (i ^ oo) , 



(A.7) 



which is the expected behaviour in the one-dimensional limit. 

For the rectangle with fixed boundary conditions one can use the same procedure, but 
in this case, owing to the fixed boundary, the Green function GrxT, is no longer transla- 
tionally invariant. For the e-dependent part of the Green function it is advantageous to 
use the following representation 



Grxt{z, z') = G{z - z') + G{z + z') - g{z - z') - g{z + z>) , 



where Q is the periodic Green function for a torus of size 2R x 2T. Denoting with 'w'j-,^{R) 
the mean squared width evaluated at the symmetry point of the Wilson loop we have 

27rawji^{R) = 



k=2 



with 



Cfc 



log i^Rt/y/lT^e 

oo 

E ^^Ip^fc [(1 + cos(2A; arctan t) - {-t'f - 1 

+XW - l^k=2 fct(2fc+2)! ^kCk 
WW IrxT'^'^^ IrxT^"^^' ^Rxt{Z: z') ) 

3(32'^ + !)(! + t^)''^^ sin(2A; + 2) arctan t - 



(A.9) 
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(1 + 9t2)'=+i sin(2A; + 2) arctan 3t - (9 + sin(2A; + 2) arctan t/3(A.10) 

X{t) = [9f{R, T) - 3f{3R, T) - 3f{R, ST) + f{3R, 3T)] /RT 

f{x,y)= ^ log(x2 + y2) + ^ arctan I + ^ arctan I . (A.ll) 

Again the infinite sums can be resummed using Eq. ([A.5|). 



The quadruple integral in Eq. ( A. 9 ) has the following representation 
27r 



jRxT JRxT 



64 v-oo (2n-l) ^ (A ^')\ 

2^n=l \+qin-l ) V^'-'-^J 

where C is the Riemann zeta function and q = exp(— vrt). 

A check of these formulas can be obtained by evaluating the t ^ co limit of Eq. (A. 9), 
giving 

lim 27ro-u;?i^(i?) = log(2i?/7re) , (A. 13) 

t— »oo 

which coincides with the result for the infinite strip obtained by directly using the Green 
function of Tab. I. 
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Figure Captions 



Fig. 1. Variation of the mean squared width of the flux tube as a function of the ratio 
T/ R calculated in the gaussian model and compared with numerical simulations 
in the 2D XY model (open circles and triangles for periodic and fixed boundary 
conditions respectively) and the 3D ^2 gauge theory (black dots). The continuous 
and dashed lines correspond to periodic and fixed boundary conditions, respectively. 
The dotted line indicates the infinite strip limit. 

Fig. 2. Density of the flux tube generated by a 30 x 30 Wilson loop at /5 = 0.7460. The 
dashed line is the fit to the gaussian distribution of Eq. (^). The continuous line 
is a fit to the function described in the text while the dotted line is the plaquette 
expectation value in the ordinary vacuum. 

Fig. 3. Squared width of the flux tube in units of sigma. The open symbols are the data 
of the ^2 gauge theory and are taken from Tab. II, while the black circles are data 



from the Ising interface taken from Ref. |15| 
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